Resistance jumps and the nature of the finite-flux normal phase in ultra-thin 

superconducting cylinders 
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Recent observations have revealed the emergence of an unusual normal phase when a magnetic 
flux threads an ultra-thin superconducting cylinder. Moreover, with increasing temperature, the 
resistance rises in a series of abrupt jumps. These phenomena are explained using a novel approach, 
which allows calculation of the resistance in the presence of amplitude and phase fluctuations of 
the superconducting order parameter, and at the same time introduces a local probe of the current 
and chemical potential. It is demonstrated that phase fluctuations lead to the sequential breakdown 
of local superconducting phase correlations, resulting in the formation of normal weak links, which 
give rise to the emergence of the normal phase is a stepwise manner. Finally, specific predictions 
are made on how the experimental observations change with the cylinder geometry. 

PACS numbers: 74.25.fc, 73.23.-b, 71.10.Fd 



Almost half a century ago. Little and Parks |I] per- 
formed one of classic experiments in superconductiv- 
ity: they demonstrated that the critical temperature of 
a cylindrical superconductor (SCR) varies periodically 
with the magnetic flux $ threading the cylinder; the pe- 
riod $0 = hc/2e reflects the charge 2e of the Cooper 
pairs. This effect is well understood within the BCS 
mean-field model 2J , as the kinetic energy of the electrons 
depends periodically on the magnetic flux. In fact, it has 
been predicted [3i and subsequently observed 0] that if 
the cylinder circumference is reduced to the same order 
as the SCR coherence length, the flux can drive the SCR 
into its normal phase even at the lowest temperature. 
Surprisingly [4|, the resistance of this low-temperature 
normal phase, only weakly dependent on temperature, is 
considerably smaller than the high temperature, normal 
state resistance of that same sample. A subsequent ex- 
periment has revealed even more intriguing results: the 
resistance increased with rising temperature in a series 
of steps , which broadened with applied magnetic flux. 
The origin of these steps, and the nature of the resis- 
tive state at finite flux have become a subject of much 
debate. Originally it was suggested that the steps 
arise due to consecutive events of phase separation in the 
vicinity of a quantum phase transition (see also [6]). It 
was later demonstrated [3 [8], however, that such a sce- 
nario is inconsistent with the system parameters, and a 
spontaneous transition to a symmetry broken order pa- 
rameter is not possible. Alternatively, a mean-field-like 
calculation [71 [H] suggested that the observation might 
be due to disorder induced fiuctuations in the coherence 
length. However, it was then claimed [S] that the large 
variation in the coherence length necessary to explain 
the data was inconsistent with other features in the ex- 
periments. Thus, the experimental observations remain 
hitherto unexplained. 

In a two-dimensional system at the critical temper- 
ature Tkt phase fluctuations of the pairing amplitude 
drive vortices and anti-vortices to unbind and proliferate 
through the system [IIJ' . The destruction of global phase 



coherence drives the loss of perfect conductance ^Ij . It is 
therefore imperative for any theory that attempts to de- 
scribe the loss of superconductivity in these ultra-thin 
cylinders to take phase fluctuations into account. In 
this letter we utilize a new formulation 12 of trans- 
port through low-dimensional, possibly disordered SCRs 
(attached to two metallic leads), in the presence of fi- 
nite temperature and magnetic field, to study a micro- 
scopic model of the Little-Parks effect. This ab initio 
tool, which employs an exact formula for the current, 
takes full account of thermal phase fiuctuations of the 
superconducting (SC) order parameter, while neglecting 
its quantum fluctuations. This method has already been 
shown [12] to reproduce, for example, the classic Little- 
Parks effect and the conductance characteristics near the 
Kosterlitz-Thouless transition. Moreover, the new tool 
introduces a local probe of the normal and SC current 
and chemical potential distribution within the sample, 
giving us the opportunity to expose and understand the 
local physical processes that drive the loss of supercon- 
ductivity and concomitant steps in the resistance. In 
the present case the intermediate region, an ultra-thin 
cylinder, sandwiched between two metallic leads, is de- 
scribed by the disordered negative-f/ Hubbard lattice 
model. The parameters that characterize the model, in 
addition to temperature T and magnetic flux $, are the 
on-site attraction U (all energies are expressed in terms 
of t, the uniform hopping integral), the typical disorder 
W - the width of the distribution of the on-site energies, 
and the average density n (see [H] for more details) . The 
coupling between the noninteracting leads and the central 
region results in a length independent contact resistance, 
which can be easily dealt with. The magnetic field is 
applied parallel to the axis of the cylinder, which has a 
finite wall thickness d and inside radius r (expressed in 
terms of the lattice constant a). 

Our strategy to study the Little-Parks effect in an 
ultra-thin cylinder is to first establish the geometry and 
key parameters required to recover the main experimen- 
tal phenomena, and then compare and contrast the obser- 
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FIG. 1: (Color online) (a) The variation of resistance with 
flux <E> and temperature T, with the normal state at "1> = 
first emerging at temperature To = O.Ilt. Rq = 17h/e^ is the 
normal state resistance. The plot is perfectly symmetric in 
magnetic flux, (b) The variation of resistance with flux for 
four different cylinder wall thicknesses d. (c) Cuts through 
the resistance surface at four different magnetic fluxes. The 
horizontal gray lines denote the emergence of the first three 
steps, (d) Cuts through the resistance surface at zero flux 
for two different length SCRs. (e) The steps in the resistance 
against step number for systems of length 60a and 240a at 
$ = 0. 



vations to experiment. Secondly, we will take advantage 
of our new capability to study the microscopic observ- 
ables to highlight the underlying mechanisms that drive 
the emergence of the normal state. Finally, we will use 
our formalism to make predictions about further observ- 
able phenomena that could verify the validity of our ap- 
proach and help pin down the nature of the transition. 

In order to determine a suitable geometry for our cylin- 
der, we first note that the resistance was not a periodic 
function of the magnetic field because of flux penetrating 
through the sample walls [U [5] . Our first numerical re- 
sults in Fig. fl|b) depict the magnetic field dependence of 
the resistance R for several wall thicknesses d. We find 
that the choice d w 0.35r (which falls within the typical 
range of the experimental devices), has a reasonable re- 
semblance to the experimental results, and the numerical 
results presented hereafter are for that choice of thick- 
ness (in all the calculations presented in this Letter, we 
set the parameters U — 1.6t, W = O.lt, n — 0.9, use 
a cylinder of circumference 11a, and length L = 60a). 
Fig. [iTa) depicts the dependence of the resistance with 



flux and temperature. In addition to the expected sup- 
pression of superconductivity by temperature, it demon- 
strates the emergence of the normal phase at low tem- 
perature around $ « $o/2, and the recurrence of the 
SC phase for $o/2 < $ < 'I'o- This surface displays the 
same qualitative features as the experimental results in 
Ref. [Ij. In Fig. [l|c) we plot cuts through the surface 
at constant threading flux. The cut taken at zero flux 
contains a substantial SC phase region (where the resis- 
tance is only due to the contact resistance to the leads) 
before the resistance increases into the normal phase in 
a series of steps. Increasing flux lowers the transition 
temperature, until at half-integer flux the system starts 
out in a normal phase whose resistance then increases 
further with rising temperature. This normal state per- 
sists to almost T = at which point the weakly coupled 
superconducting (SC) regions become phase locked and 
the resistance rapidly drops. This can also possibly be 
seen in some of the experimental results [5! . 

Having outlined the main qualitative behavior we now 
focus more closely on the emergence of the steps in the 
resistance, which are also seen in the experimental re- 
sults [5] . For the cuts at different fluxes in Fig. [TJc) , the 
steps in the resistance emerge at similar values of the 
overall resistance roughly independent of flux. Interest- 
ingly, for the longer L = 240a wire shown in Fig. Iljd) 
more steps emerge. To examine these steps in more de- 
tail. Fig. Ilje) depicts the dependence of the resistance 
step value on the step number, for cylinders of two dif- 
ferent lengths, L — 60a and L = 240a. To reliably distin- 
guish the steps we focus on only the lowest steps (these 
steps can be also positively identified from the successive 
emergence of boundaries in the phase coherence plots in 
the lower panel of Fig.[2J see below). The resistance steps 
are consistent with either an exponential (as deduced in 
the experimental paper j5j) or power law dependence. 

Having reproduced the experimental results, we now 
take advantage of our ability to produce local current 
and voltage maps, and differentiate between SC (Cooper- 
pair) and normal current [12j . to expose the microscopic 
mechanism that leads to the emergence of the normal 
phase with increasing temperature and flux. We start 
at low temperature [point (i) in Fig. flVc)] well below Tq, 
where the system first assumes a finite resistance. Here, 
as we can see in the current map, [top panel of Fig. [2ji)] 
all the current is SC. This conclusion is substantiated 
by the SC order parameter (second panel) being nearly 
constant along the wire. Defining the phase of the order 
parameter 9^ at a point x along the cylinder averaged 
over the circumference, we can see in the lower panel that 
the long-range SC correlation function (cos(6':ri — 6*3,2)) is 
constant across the system, demonstrating that the entire 
SCR is phase coherent. We also note (third panel) that 
there are no vortices in the superconductor, and there is 
no voltage drop along the sample (fourth panel), so the 
only source of resistance is the contact resistance. 

As temperature increases to T = 1.2ro, [Fig. [2l^ii) , sec- 
ond panel, corresponding to point (ii) in Fig. flTb)], the 
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FIG. 2: (Color online) The current maps (upper) at four points labeled in Fig. [l[e). The normal current is shown by cyan 
pointers and supercurrent in magenta darts, whose length and direction reflects the local current flow. The map overlays the 
spatial variation of |A] shown by dark orange. The second row of plots show the column averaged variation of |A|, the third 
row the column averaged number of vortices per site, and the fourth row the local potential. The bottom row shows the column 
averaged correlation {cos(&i:i — 9x2)}, with strong correlations in white, and weak in black. 



SC correlations slightly weaken at around x/a ~ 17 for 
this specific realization of disorder. While a mean-field 
calculation would still lead to long-range coherence in 
the sample, temperature also enhances phase fiuctuations 
across the weaker SC region. These weaker local SC cor- 
relations allow nucleation of vortices at lower energy (see 
Fig. l2jii) central panel) driving a local loss of correla- 
tion across the junction. This can be seen in the lower 
panel, where the system separates into two disconnected 
strongly correlated SC regions, giving rise to an effec- 
tive Josephson junction, and a finite voltage developing 
across the junction (fourth panel). At this temperature 
the resistance effectively jumps to a finite value. The cur- 
rent through the system is naturally reduced, and in fact 
changes its character from SC to normal and back to SC 
as it crosses the junction (top panel). The bottom panel 
of Fig. l2lii) indicates also that weaker correlations are 
developing between the two sides of x/a ~ 45, which can 
be thought of as an effective Josephson junction whose 
maximal current is still larger than the current in the sys- 
tem, and thus no voltage drops across it. As temperature 
increases further, there is little change in the resistance, 
which is basically due to the normal resistance of the 
Josephson junction at x/a ~ 17, until thermal fiuctua- 
tions lead to further nucleation of vortices, voltage devel- 
oping across a second Josephson junction at x/a ~ 45, 
and a second jump in the resistance [see Fig. |2jni), cor- 



responding to point (iii) in Fig. [He)]. Increasing the 
temperature further leads to the formation of additional 
weak links, resulting in more jumps in the resistance. 

Unlike the mechanism suggested in Refs.[71 [5], one 
does not need to invoke large variations in the coherence 
length of up to 40% along the sample to quench the SC 
order. From the data in the first panel of Fig. [2]ji) we can 
estimate the local coherence length using ^ = hvp/nA 
as varying by no more than ~ 15% across the sample. 
In this low dimensional system, small variations in the 
coherence length are sufficient to induce separation of 
the system into local SC regions, connected by normal 
weak-links, due to the increasing importance of ther- 
mal fluctuations. On the other hand, the system indeed 
separates into normal and SC regions, as suggested by 
Ref. , but not due to spontaneous phase separation near 
a critical point, but due to the interplay of disorder and 
temperature driven fluctuations. Similar behavior is ob- 
served upon keeping the temperature fixed, but increas- 
ing the magnetic fiux in Fig. ^fiv) [see also point (iv) in 
Fig.lslb)]. Here the loss of coherence is not driven by fluc- 
tuations, but instead the vortices, induced by the finite 
magnetic field, disrupt the SC state even at very low tem- 
peratures. Nevertheless, we still expect that as a func- 
tion of fiux, the resistance will rise in steps, similar to the 
steps that emerge with rising temperature. This indeed 
is borne out by the numerics, as depicted in Fig. p[b). 



(a) Diameter 



(b) Constant temperature 



^ 



1 

' . . ■ ■ ' / 

,./^" (ii) Bulk SCR/' 
-; [ L breakdown ' 


\ "■■■ 

(i) LP SCR 

breakdown 


'y> 


...^ 


1 


■-V 


/ 
/ 

/ 


Normal / 


^\,. 


1 
t 


Superconducts 




(c) Disorder 



10" 

5^ 10-1 



10- 



- ir = i.2y 






/ 


7^ 


W = 0.21/ 


1 J. 


r 




vV = 


= 0.1t , 



0.5 



0.5 



1 

T/To 



FIG. 3: (Color online) The influence of cylinder geometry on 
the resistance: (a) The SCR to normal state transition for 
different cylinder diameters for a d = cylinder (green solid) 
and d = 0.35r (orange shading). The red dotted line (i) shows 
the mean-fleld prediction for the LP superconductor-insulator 
transition, and the blue dashed line (ii) shows the mean-field 
prediction for the breakdown of superconductivity due to flux 
penetrating the walls, (b) The variation of resistance with 
flux for four different cylinder wall thicknesses r. (c) The 
zero flux superconductor-insulator transition for three levels 
of disorder. 



Furthermore, the lower panel of Fig. ^iv) shows that for 
the same system, the breakdown of the SC phase cor- 
relations takes place at the same location, x ~ 17a as it 
did with increasing temperature in Fig.^iv), so pointing 
toward the same disorder driven transition. The observa- 
tion that the resistance jumps are almost identical for a 
longer sample [see Fig. IlTd)], except that there are more 
of them, further supports the scenario where each step is 
due to the addition of a single normal weak link. 

So far we have concentrated on choosing parameters 
to best fit the experimental observations. The formal- 
ism allows us, naturally, to explore any range of param- 
eters, and make specific predictions. We have already 
reported in Fig. [TJb) the flux dependence of the resis- 
tance for cylinders of different thicknesses. As expected, 
the thicker the cylinder walls, the larger the magnetic 
field penetrating the actual sample, and the less peri- 
odic the signal is, until &t d ^ 1.05r, the SC phase 



does not reemerge beyond $ = <I>o/2 (the actual value 
of this thickness may depend on disorder) . To study this 
in more detail Fig. IsFa) depicts the full phase depen- 
dence on wall thickness. For a two-dimensional cylinder 
(d — 0) the numerical separatrix between the SC and 
normal phase deviates from the mean-field expectation 
r = ^0 niin„ |n-|-<i>/$ol [labeled (i)] at small cylinder radii 
r, due to the fact that in such quasi one-dimensional sys- 
tems, thermal fluctuations are sufficient to disrupt the SC 
phase even at zero flux. The filled orange area depicts 
the SC phase for a finite-thickness cylinder d = 0.35r. At 
small flux we see the same periodic variation as the d = 
cylinder. However, flux penetrating the walls will itself 
destroy the SC phase along [3] r = I^^^J^q^JX-T jT^ 
[labeled (ii)], with ^q the SC coherence length used as 
a fitting parameter, and so place a bound on the phase 
boundary. 

We have also shown in Fig. [iTd) that the values of 
the resistance plateaus depend only weakly on the cylin- 
der length, though a longer cylinder is predicted to ex- 
hibit more steps. Additionally, Fig. IsFb) demonstrates 
that one should also see steps in the resistance as a func- 
tion of flux, for a flxed temperature. Interestingly, these 
steps persist even above the critical temperature, which 
is consistent with the picture presented above - the criti- 
cal temperature corresponds to the formation of the flrst 
normal weak link, but additional normal links are formed 
with increasing flux. Finally, in Fig. Islc) we plot the 
variation of resistance for two higher levels of disorder. 
At higher disorder amplitude the transition temperature 
decreases, and the normal state resistance increases. In- 
terestingly, the resistance steps persist even for a highly 
disordered cylinder {W = 1.2t) that exhibits a normal 
resistance at low temperature. This prediction could be 
readily checked experimentally. 

To conclude, we have used ah initio simulations to ex- 
amine recent experiments on the quantum Little Parks 
effect. The simulations demonstrate step- wise destruc- 
tion of the SC phase, stemming from phase fluctuations 
breaking down SC coherence in those parts of the cylin- 
der that have a weaker BCS order parameter, due to 
disorder. The formalism also allowed specific predictions 
of further phenomena that could verify this hypothesis. 
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